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1. Introduction and Background 

1.1. Some Basic Issues in Relativistic Quantum Information 

We assume the readers are somewhat familiar with the Unruh effect [1] on the one hand 
and the basic issues of quantum information on the other [2] and will only highlight 
the relativistic aspects of both of these topics here. First, the quantum field acting as 
an environment to the discrete (qubits) or continuous variables (oscillators) in quantum 
information processing - we will refer to these point-like physical objects with internal 
degrees of freedom as detectors (Unruh-DeWitt detector being the familiar one [3]) - 
will necessarily exert environmental influences on the system. Second, the motional 
states of the detectors (e.g., inertial or accelerated, uniformly or otherwise) will affect 
both the quantum decoherence and entanglement dynamics of these detectors. 

1.1.1. Quantum Field Effects The presence of a quantum field is unavoidable, as it 
acts as an ubiquitous environment to the qubits or detectors in question. Two basic 
issues of quantum information need to be included in one's consideration are: 
Quantum decoherence Coupling to a quantum field can induce decoherence of a 
single qubit or oscillator, but their mutual influences mediated by a field can lessen the 
degree of decoherence if the two qubits are placed in close range [I]; 
Entanglement dynamics The entanglement between two qubits or oscillators 
changes in time as their reduced state (after coarse-graining over the field) evolves; 
it also depends on their spatial separation [U [5] . 

1.1.2. Kinematical effects 

Unruh effect A uniformly accelerated detector coupled with a quantum field in the 
Minkowski vacuum would experience a thermal bath of the field quanta at the Unruh 
temperature proportional to its proper acceleration. This was first discovered by Unruh 
using time-dependent perturbation theory pQ. Generalized considerations follow in the 
works of Higuchi et al [6] and Louko et al [TJ. Exact solutions going beyond these test- 
field descriptions were found by Lin and Hu [HI E] with interesting new physics. 
Non-uniform acceleration The kinematical viewpoint has proven to be more 
malleable and adaptable than the traditional geometrical (global concepts like event 
horizon) viewpoint. We will mention how newer models in the 90s such as the RSG 
model [10], especially those which incorporate open quantum system concepts such as 
the RHAK models [TT1 IT2] . have aided in treating non-uniform acceleration, in work 
from the 90s (e.g., [13]) to now [T3] . 

Mutual influences The influence of one detector on the field will propagate in space 
and affect other detectors after some time. These causal mutual influences propagating 
back and forth is a source of non-Markovianity in multi-detector theories. They act to 
augment the quantum coherence between two detectors placed in close range. Another 
source of non-Markovianity is the long range autocorrelation of the quantum field. 
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1.1.3. Relativistic Effects Furthermore, objects in a relativistic system may behave 
differently when observed in different reference frames, so we have: 
Frame dependence Quantum entanglement of two objects localized at different 
positions on a spatial hypersurface is a kind of spacelike correlation, the time evolution 
of their entanglement will depend on how the spacetime is foliated by spacelike 
hyper surf aces. 

Time dilation For moving objects with worldlines parameterized by their proper 
times, their time dilations observed in a reference frame will naturally enter the dynamics 
in that frame. 

Projective measurement local in space Quantum states make sense only in a given 
frame where a Hamiltonian is well defined [15J. Two quantum states of the same system 
with quantum fields in different frames are comparable only on those totally overlapping 
time-slices associated with certain moments in each frame. By a measurement local in 
space, e.g. on a point-like UD detector coupled with a quantum field, quantum states 
of the combined system in different frames can be interpreted as if they collapsed on 
different time-slices passing through the same measurement event. Nevertheless, the 
post-measurement states will evolve to the same state up to a coordinate transformation 
when they are compared at some time-slice in the future. In a two-detector system with 
the first detector being measured at some moment, the reduced state of the second 
detector collapsed in different frames will become consistent once it enters the future 
lightcone of the measurement event. [U [161 El HE]- 

1.2. Unruh effect via perturbation theory and quantum information via exact solutions 

Time-dependent perturbation theory (TDPT) was used by Unruh originally to show the 
detector response to uniform acceleration. For a comprehensive description of Unruh 
effect, see, e.g., the recent review of [19]. With the infusion of quantum open system ideas 
in the 90s these TDPT results were later found to hold only in the Markovian regime, 
corresponding to the limits of ultra-high acceleration or ultra-weak coupling. Discovery 
of exact solutions in the 2000s showed that the transition probability calculated from 
the infinite-time TDPT is valid only in transient under restricted conditions. We will 
develop the perturbative theory further in Section 2 and comment on these developments 
in the last two sections. In Section 3 we will introduce two other more general models for 
moving detector- quantum field interaction, namely, the RSG and the RHAK models, 
treating the detectors as harmonic oscillators rather than the two level system as in the 
original Unruh derivation. We will also bring in the broader scope provided by the theory 
of open quantum systems exemplified by the quantum Brownian model, where the use 
of reduced density matrix and influence functionals opens the way to exploring the full 
parameter range of detector-field interaction including self-consistent backreaction. This 
opens the door for quantum information inquires. In Section 4 we describe the detector- 
field dynamics from the exact solutions of one such model, and identify the limitations of 
TDPT. In Section 5 we give an important example of relativistic quantum information, 
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that of quantum teleportation which uses pretty much all of the relativistic and quantum 
information elements developed, such as frame dependence and entanglement dynamics. 
We end with some suggestions on further developments. 

2. Nonstationary detector within first-order perturbation theory 

In this section we summarise recent results about the transition rate of a pointlike 
detector within linear perturbation theory, in situations where neither the detector 
trajectory nor the state of the quantum field is assumed stationary [TJ EH EQ 122] . 
The central issue is to isolate switch-on and switch-off effects from phenomena that are 
genuinely due to the acceleration and to the state of the field. 

2.1. Transition probability 

We consider a pointlike two-level detector that moves in a spacetime of dimension d > 2 
along the worldline z(r), where the parameter r is the detector's proper time. The 
motion is driven by an external agent who is decoupled from the detector's internal 
degrees of freedom and from the quantum field to which the detector couples. 

The detector's internal Hilbert space here is two-dimensional, spanned by the 
orthonormal basis states |0)<f and \oj)d whose respective energy eigenvalues are and u, 
with oj 7^ 0. For u > 0, |0)d is the ground state and \cj)d is the excited state; for 
u> < 0, the roles of the states are reversed. A generalisation to a countable number of 
nondegenerate energy eigenstates would be straightforward. 

The spacetime contains a free real scalar field <p, whose mass and curvature coupling 
parameter may be arbitrary. The detector is coupled to <fi linearly, by the interaction 
picture Hamiltonian 



where A is the coupling constant and Q is the detector's monopole moment operator. 
The switching function \ specifies how the interaction is turned on and off. We assume 
X to be smooth, nonnegative and of compact support. We also assume the trajectory 
z(t) to be smooth. 

We denote the initial state of the field by |-0 O )> an d we assume |^ ) to be regular in 
the sense of the Hadamard property [231 EI] • The detector is initially prepared in the 
state \Q)d- 

We work within first-order perturbation theory in A. After the interaction has 
ceased, the probability for the detector to be found in the state \u)d, regardless the final 
state of the field, is [251 126] 



int 



Ax(r)Q(r)0(z(r)) , 
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and the correlation function W(t',t") := (^|0( z ( r/ ))0( z ( r "))l' ? / ; o) is the pull-back of 
the Wightman function to the detector's worldline. The prefactor A 2 |d(0|Q(0)|c<;)d| 2 
in (EJ) depends only on the detector's internal structure, while T (cu) ([3]) encodes the 
dependence on |^o)) the trajectory and the switching. With minor abuse of terminology, 
we refer to T (u) as the transition probability. 

2.2. Transition probability without distributional integrals 

While formula ([3]) for the transition probability is as such well defined, it is not well 
suited for discussing how the probability depends on the switching function, especially 
when the switching becomes sharp. The correlation function W is not a genuine function 
but a distribution. When W is represented by a family W e of functions that converge to 
W as e — > 0+, the sense of convergence entails that the limit e — > + is taken in ([3]) only 
after the integrals are evaluated [271 [28| 1291 [30] . The sharp switching limit may hence 
not necessarily be brought under the integrals and the e — > 0+ limit in (jHJ) [3U [32], 133] . 

What is needed is to re-express ([3]) in terms of the genuine function Wq : = 
lim e _i.o + W e , where the limit is understood pointwise. The results for d = 2, d = 3 
and d = 4 are [7J US I2Q 122] 



F d=2 {u) =2 Re / du X (u) ds X (u - s) e~ luls W (u, u - s) , (4) 

1 poo poo POO 

T d=3 (u) =- du [ X (u)] 2 + 2 / du X (u) / ds X (u - s) Re [e'^Woiu, 

^ J —oo J —oo Jo 



U — S) 



/oo 2 poo 

dw[x(n)] 2 + — / — / d«x(«)[x(«)-x(«-*)] 



+ 2 y dw ^ ds X ( M - s) Re ^ s W (u, u - s) + _) , (6) 

and those for d — 5 and c? = 6 can be found in [22] in the special case of a Minkowski 
space massless field in the Minkowski vacuum. The crucial point is that in addition 
to an expected integral term that involves Wo, there are also additional terms that 
depend on the switching. These additional terms are remnants of the distributional 
singularity of W, and they are absent only for d = 2, where the singularity of W is 
merely logarithmic. 

The Hadamard property of the Wightman function guarantees that the integrals 
in (SJ) — (J6j) are convergent at s = 0. We assume that any singularities of Wo at s > 
are integrable. Such singularities can occur for example when the spacetime has spatial 
periodicity so that points on the detector's trajectory can be joined by null geodesies 
that circumnavigate the space [22J. 
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2.3. Transition rate 

When both the detector trajectory and the quantum state of the field are stationary, 
in the sense that they are invariant under a Killing vector that is timelike in a 
neighbourhood of the trajectory, a transition rate per unit time may be defined by 
making the switching function time-independent and formally factoring out the infinite 
total time of detection [H El HH [251 I2S1 [Ml EH EH EZl EH E9] . In time- dependent 
situations this procedure is however not available, and separating the switching effects 
from the acceleration effects becomes delicate p[I3E3llJlimil21l33Sll35]. 

To define a transition rate in the nonstationary setting, we consider the limit in 
which the detector is switched on an off sharply. We let the switching function x take 
the value unity from proper time To to proper time r, where r < r, and we assume 
that the switch-on takes place over an interval of duration 5 before To and the switch-off 
takes place over an interval of duration 5 after r, in a manner discussed in [7J [20]. The 
limit of sharp switching is 5 — > 0. 

We regard the response function J 7 as a function of the switch-off moment r, and we 
define T T := dF /dr. T T may be regarded as the detector's instantaneous transition rate 
per unit proper time, observationally meaningful in terms of consequent measurements 
in identical ensembles of detectors [TJ. 

For d = 2 and d — 3, taking the the 8 — > limit in (JJJ) and fl5]) is immediate and 
yields a finite result for the transition probability. For d — 4, the 5 — > limit in ([6]) 
contains a divergent term proportional to In § [7J [20]. This divergent term depends on 
the details of the switching but it is constant in time, and it is also independent of the 
trajectory and of the quantum state. The divergent term does hence not contribute to 
the transition rate. Physically, the 5-^-0 limit means that we take the switching to be 
rapid compared with the overall duration of the interaction: focusing on the transition 
rate allows us to discard from the transition probability the numerically dominant piece 
that only depends on the details of the switching. Collecting, the 5—^0 transition rates 
for d = 2, and d = 3 and d = 4 are given by 

d = 2: F T (u) = 2Re ds e~ luJS W (t, t - s) , (7) 

Jo 

d = 3: P T (u) = - + 2 f ds Re [e~ iuJS W {T, r - s)] , (8) 
4 Jo 

d = 4: ^( w ) = _^ + 2 | aT dsRe (e-W„(r,r- S ) + ^) + ^ , (9) 

where At := t — t . (J7J) and (jHJ) are valid as 5 — > at fixed A, provided A is so small that 
the total transition probability remains within the validity domain of the perturbative 
treatment. ([9]) is valid as 5 — > provided A simultaneously approaches zero so fast 
that it is bounded in absolute value by k/ y | In 5\, where the positive constant k is so 
small that the total transition probability remains within the validity domain of the 
perturbative treatment. 
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For d = 5 and d = 6, we specialise to a massless field in Minkowski spacetime in the 
Minkowski vacuum [21] . The transition probability contains again a term that diverges 
as S — > 0. For d = 5 the divergent term is constant in time, and the transition rate has 
the finite 5 — > limit 

, _ ^ , N 4o; 2 + z 2 (V) 1 /" Ar , / sin(ws) w \ w 
a = 5 : y- r fa;) = - 




647T 47T 2 ./n I /r /a x2l3 S 2 / 47T 2 Ar ' 

(10) 

where Az := z(u) — z(u — s). For d — 6, by contrast, even the transition rate contains 
a term that diverges for generic trajectories as 5 — > 0, proportionally to z • z\n5. 
This means that the divergences due to the rapid switching cannot be isolated from 
the acceleration effects for d = 6. The sole exception occurs for trajectories whose 
scalar proper acceleration a/^ 2 is a constant, including as a special case all stationary 
trajectories. For such trajectories the d = 6 transition rate remains finite as 5 — > and 
is given by 

, R oo{oo 2 + -± 2 ) 1 [ At , /cosM 1 . 3o; 2 + z 2 
a = 6: Jvfa;) = — + — - / as 



24vr 2 2n 3 J V[(Az) 2 ] 2 s 4 6s 2 

3a; 2 + z 2 1 
12tt 3 At 6tt 3 At 3 ' 



(11) 



S.^. Applications 



When both the detector trajectory and the quantum state of the field are stationary, the 
transition rate formulas (ITjl-f lllj) reduce to the well-known formulas in which stationarity 
is assumed at the outset [25], [26], [38] - We re-emphasise, however, that formulas (J7j)- (fTTj) 
apply in genuinely time-dependent situations. 

A showcase example is a Minkowski spacetime trajectory that is asymptotically 
inertial at early times and of asymptotically uniform linear acceleration at late times, 
with the field in the Minkowski vacuum. Within the perturbative treatment, the 
transition rate is duly found [33] to interpolate between that in inertial motion and 
that in uniform linear acceleration, describing thus the onset of the Unruh effect pQ. 

Other applications can be found in [7J [201 1221 [331 US]- Slowly- varying acceleration 
is discussed in [471 HH] ■ 



2. 5. Other definitions of the transition rate 

To end this section, we mention two alternative definitions of the transition rate. 

First, the transition rate of a pointlike detector in flat spacetime can be defined by 
first giving the detector a spatial size, specified covariantly in terms of the detector's 
instantaneous rest frame, and at the end taking the pointlike limit [3T1 1321 [38] . The 
results agree with those obtained via smooth switching in the common domain of 
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validity [2TJ EEJ E2J ESI SHI ED]- A related procedure that replaces spatial size by a 
pole prescription in proper time is discussed in [37]. 

Second, a spatially extended detector in flat spacetime can be reinterpreted as a 
pointlike detector with an energy cutoff that is specified covariantly in the detector's 
instantaneous rest frame j3HJEQ]. Definition of the transition rate via this energy cutoff 
can be generalised to curved spacetimes at least when the spacetime has a sufficient 
amount of symmetry [33 |50] . 

3. Detector-Field Interaction: UD, RSG, RHAK models 

Quantum mechanics for the single particles in some non-linear potentials such as 
anharmonic oscillators, Morse potential, etc. are exactly solvable. But in field theory, 
since a field has infinitely many degrees of freedom, a small non-linearity can create huge 
difficulty in calculations. One can at most do perturbation theory or self-consistent 
approximations around some non-trivial background field configuration, where the 
calculation involves essentially Gaussian integrals. However, if the potential of the 
detector is that of a harmonic oscillator (HO) and the quantum state is in a Gaussian 
form, it is possible to solve the full dynamics of the combined system of the detectors 
and the field non-perturbatively. 

3.1. Raine-Sciama- Grove (RSG) Model 

A notable non-perturbative model is that of Raine, Sciama, and Grove pU], where they 
found the late-time expectation values for the stress tensor of a massless scalar field in 
(1+1 )D Minkowski space. The method is generalized to the case with a point-like UD 
detector in a massless scalar field in (3+1 )D Minkowski space, and the whole history of 
the combined system is solved in [8]. 

3.2. Proxy to Quantum Brownian Motion Models 

Suppose the internal degrees of freedom of the UD detector are HOs. Then the combined 
system of N UD detectors and a quantum field is an (N + oo)-HO system, which 
is linear and exactly solvable. Unruh and Zurek [51] have studied a model where 
a hamonic oscillator interacts with a massless scalar field in 2-D. They derived the 
exact master equation for the reduced density matrix of the system (oscillator) at a 
temperature determined by the initial state of the field, and observed some general 
features different from the conventional Markovian results valid for an ohmic bath at 
ultra-high temperature made known earlier in the famous paper of Caldeira and Leggett 
[52] . One feature is the dependence of the ultraviolet cut-off in the master equation and 
the reduced density matrix, and thus also in the von Neumann entropy of the system. 

More general non-Markovian behavior was explored by Hu, Paz and Zhang [53] 
who derived an exact master equation with nonlocal dissipation and colored noise for 
the system of one harmonic oscillator (detector) interacting with a thermal bath of 
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n-harmonic oscillators. For quantum decoherence they identified the low temperature, 
supra-ohmic regime as a noticeable departure from the Markovian behavior (see followup 
in [M])- Using this model as a theoretical tool with the help of quantum open system 
ideas, many of the basic issues we listed in the beginning can be addressed effectively. 

Generalizing a bath of n-harmonic oscillators of time- dependent frequencies to a 
quantum field was subsequently done by Hu and Matacz [55] for moving detectors- 
quantum field interactions. Because the treatment is given in quantum optics language 
- the quantum states described by the squeeze, rotation and displacement operators 
and the dynamics in terms of parametric amplification - their results are immediately 
applicable to "atomic- mirror" -optical systems [56]. It also served the intended purpose 
of bringing open systems methods and concepts to quantum field theory. The influence 
functional treatment they used incorporates the backreaction of the environment on the 
system (which could be either the quantum field or the harmonic oscillator depending 
on what one is after) in a self-consistent way. In particular they showed how the Unruh 
and Hawking temperatures can be identified from the noise kernel using this method. 

Viewing the Unruh effect from this perspective, since in the QBM model there are 
nontrivial activities at zero temperature j5TJ |53j |55], we note that even for the zero 
acceleration a = case the detector is not just laying idle but has interesting physical 
features due to its interaction with the vacuum fluctuations in the quantum field. 

3.3. Raval-Hu-Anglin-Koks (RHAK) Models 

A model of N detectors in arbitrary relativistic motion interacting with a common 
quantum field (but not with each other) was proposed by Raval, Hu, Anglin [TT] . They 
calculated the influence of quantum fields on the detectors in motion, and the mutual 
influence of detectors by the action of fields via the Langevin equations derived from 
the influence functional. They introduced the notion of self and mutual impedance, 
advanced and retarded noise, and the new relations between noise-correlations and 
dissipation-propagation. They show the existence of general fluctuation-dissipation 
relations, and for trajectories without event horizons, correlation-propagation relations. 
Raval, Hu and Koks [12] used this model to explore different trajectories of the moving 
detectors in a quantum field and showed that this is a more feasible way (over the 
traditional global geometric view which relies on the existence of event horizons) to 
address situations where the spacetime possesses an event horizon only asymptotically, 
or none at all. Examples studied there include detectors moving at uniform acceleration 
only asymptotically or for a finite time, a moving mirror, and a two-dimensional 
collapsing mass. They show that in such systems radiance indeed is observed, albeit not 
in a precise Planckian spectrum. The setups in this model have been adopted in the 
study of charge particle motion by Johnson et al [57] in an electromagnetic field and by 
Galley et al [58] for the self-force of masses moving in a gravitational field. 
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3.4- Moving Detectors- Quantum Field Interaction 

Consider a model with N identical point-like Unruh-DeWitt detectors with the internal 
degrees of freedom represented by harmonic oscillators with mass mo and natural 
frequency Q, moving in a quantum field in (3+l)D Minkowski space. Here we follow 
the treatment in [16]. The action of the combined system is given by 



where a = 0,1,2,3, g aa i = diag(— 1, 1, 1, 1), d = A,B,C, - ■ ■ denotes the names of 
the detectors, d d = d/dr d , r d is the proper time for detector Qd and z d (Yd) is the 
trajectory of detector d. The scalar field $ is assumed to be massless, and A is the 
coupling constant. We consider a massless scalar field here because it is simpler and a 
good representation of the electromagnetic field. In fact all kinds of fields, massless or 
massive, bosonic or fermionic, can be considered, depending on the physics one aims at. 
The detectors do not have to be uniformly accelerated (e.g., [12]) or at rest. However, 
the motion of the detector here are assumed to be controlled by external agents, in other 
words, the trajectories or worldlines of the detectors are prescribed and not dynamical. 
If the motion of the detector becomes dynamical, it is extremely hard to get analytical 
results even in classical theory (for example, a relativistic charge in a field in classical 
[5D] and quantum field theory [OD]), since including backreaction of the field on the 
detector will alter its (test-field) prescribed trajectory. Trajectories of charged particles 
[57] and even extended objects [61] determined by their interplay with the quantum field 
have also been studied before using the influence functional method which is particularly 
suited to treating consistent backreaction effects. 

The first noticeable attractive feature of (Tl2l) is that this model is linear and thus 
easy to treat. It is arguably the simplest model for an "atom" -field interacting system 
but complex enough to give nontrivial results and insights. By "atom" here we refer 
to a spatially-localized physical object with internal degrees of freedom. By field, we 
categorically refer to dynamical variables which can be non-local in space In some 
simple setups analytic results can be obtained in the whole parameter range, and back- 
reaction to both the atom and the field can be fully studied with the help of quantum 
open systems techniques. 

Moreover, in a relativistic setting, such as for uniformly accelerated detectors or 
black holes, event horizons for the detectors can be sharply defined since the detectors 
are always localized. Also since the detectors are pointlike, they are allowed to 
be parametrized by their own proper times, which are invariants under coordinate 
transformations. This greatly simplifies the calculations in different reference frames 
when the related physics correspond only to the two-point correlators of the detectors 



S 




(12) 



| 'Nonlocal' is in the sense used by the atomic-optical quantum information community. Of course 
quantum field theory is local. 
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parametrized by two proper times. Note that this is even plausible for extended objects 
in the spirit of effective field theory [6T| [62] . 

For a uniformly accelerated UD detector in (3+l)D with proper acceleration a the 
Unruh effect [H [3j [191 125] attests that it should behave the same way as an inertial UD 
detector in contact with a thermal bath at Unruh temperature Tjj, or more precisely, 
as an inertial harmonic oscillator in contact with an Ohmic bath at Tjj [63]. However, 
examining this from the vantage point of the exact solutions we obtained, we see the 
above statements are accurate only at the initial moment. After the coupling is switched 
on, the quantum state of the field will have been changed by the detector, so the field 
is no longer in the Minkowski vacuum and it does not make exact sense to say that the 
detector is immersed in a thermal state (or any state defined in the test-field description, 
i.e., where the field is assumed not to be modified by the presence of the detector). 

A theorem by Bisognano and Wichmann (BW) [61] states that the Minkowski 
vacuum, which is uniquely characterized by its invariance under all Poincare 
translations, is a Kubo-Martin-Schwinger (KMS) state with respect to all observables 
confined to a Rindler wedge. It does not apply here because the BW theorem refers to 
the vacuum state of a quantum field alone, not the combined detector-field system. Even 
when the combined system is in a steady state, the quantum state of the interacting field 
is not invariant under spatial translations in Minkowski space, hence does not subscribe 
to the assumption of the BW theorem pertaining to Poincare invariance. Actually the 
Planck factor in, for example, (Eq.(60) in Ref.[8]), 

is a consequence of the BW theorem. Nevertheless, it is derived from only the free-field- 
solution part of the complete interacting field. Here the factor is not distorted by the 
interaction simply because the field is linear and the coupling is bilinear. For nonlinear 
fields or couplings it would have a nonPlanckian spectrum and the departure from the 
conventional picture would be more pronounced. 



4. Nonperturbative Detector-Field Dynamics 



Nonperturbatively solvable models such as ( fl2l) are particularly useful for examining 
the full features of a system which perturbative theories miss or misrepresent. They 
are essential for understanding new physics such as that associated with quantum 
entanglement whose dynamical behavior we don't really have a complete or accurate 
knowledge about. We now continue to develop the model ffT2]) . 

A quantum state of the combined detector- field system can be described by the 
density matrix p[(Q, $ x ), (Q', Q'^); x°] or equivalently, the Wigner function [511 [53] 



H/[P,E;x c 



A A n 



(14) 



§ Here we write (Q, $ x ) = £ - (A/2) and (Q', = £ + (A/2) with the boldface letters £ and A 
denoting the vectors in the configuration space. J 2? A and J 2?£ are functional integrals. 
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If we start with a Gaussian state, by virtue of the linearity of the combined system 
(TT2]) . the quantum state will always evolve in a Gaussian form in its entire history. 
Thus solving the dynamical equations for the Wigner function boils down to solving the 
time-dependent factors in the Wigner function. 

Since the field variables at some moment x° are defined on the whole time-slice 
associated with x°, the density matrix or Wigner functions at x° is also defined on 
that time-slice. In the Schrodinger picture, the evolution of a density matrix (Wigner 
function) is governed by the master equation (Fokker-Planck equation). It is possible to 
solve these equations for Gaussian states directly in simple cases (e.g. [651 [66]). However, 
when the degrees of freedom of the density matrix is large or even infinite, it becomes 
very difficult to solve the coupled equations, as the dynamics is often non-Markovian 
(the master equations for the reduced state) or non-linear in appearance (the dynamical 
equations for the time-dependent factors in the Wigner functions). Moreover, even the 
solutions are obtained and can be expressed formally, the factors in the Wigner function 
are inverse matrices with infinite dimension, which are computationally challenging. 

To get rid of these difficulties it is convenient to apply the (K, A)-representation [51] 
(or called the Wigner characteristic function [67]), which is a double- Fourier-transformed 
function of (thus equivalent to) the usual Wigner function, 

A „ A 



p[K,A;x°] = 

-<fi„>A" 



E - 



x 







= exp 


-( 

h v 



1 

2ff 



A^A 1 



(15) 



where we denote Qd and by <3> d and IT d , respectively (Pd, II X are conjugate momenta 
to Qd, $ x ), [i, v — {d} U {x} run over all the detector- and field- degrees of freedom 
defined on the whole time-slice, and the time-dependent factors Q^{x°), V j 



I'" 



x 



and 

T^fj,u(x°) are exactly the symmetrized two-point correlators (A,B) = (AB + BA)/2 of 
the dynamical variables evaluated on the x°-slice, for they are obtained by, e.g., 

ihS h5 



5A» %bK v 

= n 



p[K,A;x L 



K 



(16) 



A=K=0 



(li H ). Note that in ffT5"]) , Q^ v and V^ v are 



where <^ = $^ - and SU^ - 

defined as symmetric matrices, but IZ^ ^ 7Z UI1 in general. 

The reduced states of the system with the environment integrated out are simple 
for Gaussian states in the (K, A)-representation since a Gaussian integral gives another 
Gaussian function. For example, the reduced state of detector A with the field and 
other detectors integrated out reads 



p R [K 



A A A :x°] 



exp 



2h 2 



~ h ( {Qa)K A - (P A )A A 



{K A Q AA K A - 2A A K AA K A + A A V AA A A ) 



.(17) 



Thus, looking at the evolution of the Gaussian state (1T5T) or (fTTj) is equivalent to 
looking at the dynamics of those symmetrized two-point correlators, which would be 
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obtained more easily in the Heisenberg picture. 

4-1- Correlator dynamics for factorizable initial states 

For mathematical convenience (and to a large extent reflective of not uncommon physical 
situations) one often assumes that the initial state at Xq = to in the Minkowski frame 
is a product state of the Minkowski vacuum of the field (which is Gaussian) and the 
Gaussian state of the detectors A, B, ■ ■ ■. The detector part can be a product of the 
ground states and/or single-mode squeezed states, a multi-mode squeezed state, or any 
mixed state in the Gaussian form. The field part can be easily generated to a thermal 
state, which is Gaussian, too. 

By virtue of linearity in ( TT2l) . the operators of the detectors and the field in the 
Heisenberg picture will evolve to a linear combination of all the detector operators Qd ; Pd 
and the field operators $k ; ftk defined at the initial moment t - Then each symmetrized 
two-point correlator of the detectors for the factorizable initial state p $x (g> p d will split 
into a sum of the a-part and the v-part [8]. The a-part corresponds to the initial state 
of the detectors, while the v-part corresponds to the response to the field vacuum |0 M ). 

4-2. Divergences 

There are two sources of the divergences for the correlators in this detector-field model. 

1. In this linear system the mode functions satisfy the classical equations of motion 
(the only difference is the initial conditions). Thus they suffer the same divergences as 
the classical ones: In the UD detector theory in (3+l)D, the retarded field sourced by a 
pointlike detector diverges right at the position of the detector. One needs to introduce a 
cutoff A to regularize the 5-function in the interaction Hamiltonian between the detector 
and the field, expand the relevant mode functions of the detectors in series of A, absorb 
the divergent terms by some parameter of the model (in our model it is the natural 
frequency of the detectors Q [8], in other model it could be the mass of the detector m 
or the coupling constant A), then take the A — > oo limit to eliminate the 0(A _1 ) terms. 
The 0(A°) terms will survive after taking the limit; it gives the radiation reaction such 
as the dissipation term (~ 7Q) as those in our model or the higher derivative term as 
the one in the Abraham-Lorentz-Dirac equation [59] [JJ. 

2. The second kind of divergences are the UV or IR divergences arised in the mode sum, 
or equivalently, the divergences arised in the coincidence limit of the Green's functions 
(UV), or in the integration over the whole position space (IR). The regularization of 
the UV divergences should be consistent with those for the first kind. This kind of 
divergences is also a common feature for quantum Brownian motion. 

|| Note that there is no such divergence in the RSG model in (1+1)D [10], where the retarded field is 
regular everywhere and simply offers the radiation reaction in the equation of motion of the detector. 
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4-3. Range of validity of perturbative results 

For a single Unruh-DeWitt (UD) detector moving in (3+1) dimensional Minkowski 
space, The total action is given by ( |T2l) with d = A. Let us denote Q = Qa- Suppose 
the initial state of the system at r is a direct product of the ground state for Q and 
the Minkowski vacuum for $. It is straightforward to write down the reduced state of 
the detector in the (K, A)-representation. However, to compare with the perturbative 
results, we look at the reduced density matrix p R (Q, Q'; r) directly Transformed to the 
representation in the basis of energy eigenstates for the free harmonic oscillator Q, the 
transition probability from the initial ground state to the first excited state then reads 



o R 



h[{P 2 ){Q 2 )-(P,Q) 2 -(tf/4)} 



{[(P 2 > + (hm n r /2)} [(Q 2 ) + (h/2m Q r )} - (P,Q> 2 } 3/2 
where Q r is the renormalized frequency [8]. Expanding the symmetrized two-point 
correlators of the detector in terms of the coupling strength 7 = A 2 /(87rm ), the 
approximate value up to the first order of 7 becomes 



R I r /> a 1 ^ 



77 Ai + A - 21n(a/fi r 



(19) 



when rj = r — tq ^> a . Here Ao and Ai are large constants introduced by 
regularization (the second kind in Section 14. 2p . We see that the first term of (|19p 
gives the conventional transition probability from TDPT over infinite time. Only when 
Q r rj 3> Ai, Ao, or a is extremely large, can the second term in ( 1191) be neglected. Hence 
the conventional transition probability or transition rate is valid only in the limits of 
(a) ultra-high acceleration (a ^> Q r and Ai <C at] <C cry" 1 ) or (b) ultra- weak coupling 
(a _1 ,fl~ 1 Ai <C rj <C 7" 1 ). Only in these limits the effective temperature obtained by 
diagonalizing p R %n is very close to the Unruh temperature Ty = ha/2nkB [9], and the 
thermal bath is only slightly affected by the back reaction from the detector to the field, 
namely, the detector acts essentially as a test particle in the field. 

Note that, in obtaining ( I19j) . we have assumed a -1 <§C 77 7 -1 , when the system 
is still in transient. Indeed, in figure [1] TDPT works well only in the middle plot (with 
O(10 _1 ) < 77 < O(10~ 2 /7) for 7 = 10~ 6 and a = 6). If a < 7, the conventional transition 
probability has no chance to dominate at all. In particular, the a = case is beyond 
the reach of TDPT over infinite time, and the conventional wisdom from perturbation 
theory that no transition occurs in an inertial detector is untenable. In contrast, the 
evolution of pf \ with a = behaves qualitatively similar to those cases with nonzero 
acceleration [8]. This agrees with the expectation from the observation that the UD 
detector theory is a special case of quantum Brownian motion [511 [531 EH] • 

Ao corresponds to the time scale of switching on the interaction, namely, the "In 5" 
in d = 4 case in Section 12. 3^ so it could be finite in real processes. This implies that the 
A terms in all two-point functions will be damped out so that A will not be present 
in the late-time results. Actually (fl9l) is formally identical to the first-order transition 
probability from TDPT for a UAD with finite duration of interaction and A and Ai are 
formally the same as the divergences found in [13]. In [6] it has been shown that these 
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Figure 1. Comparison of the exact pf 1 (solid curve, Eq. (fT8]) ) and the perturbative 
result (dashed line, the first term in (|19[) ) in different time scales of 77 = r — To 
for a uniformly accelerated detector initially in its ground state and coupled to the 
Minkowski vacuum of the massless scalar field after To. The parameters in these plots 
are 7 = 1CT 6 , a = 6, il = ^/fl* - 7 2 = 2.3, A = Ai = 20, m = 1, and h = c = 1. 

divergences can be tamed if one switches on and off the interaction smoothly, so can 
A . Nevertheless, here we are looking at the real-time causal evolution problem ("in- 
in" formulation) rather than a scattering transition amplitude ("in-out" formulation) 
problem, and we never turn off the coupling, so Ai is a non-zero constant of time. 

Ai should not be absorbed by any parameter or subtracted from any physical 
quantity of this theory for more reasons: (a) The UD detector theory is not a 
fundamental theory to meet the renormalizability requirement, and the presence of 
cut-offs as physical parameters is an expected feature which characterizes the range of 
validity of this effective theory, just like the Compton wavelength of the electron serving 
as a cut-off in quantum optics; (b) Ai is not present in the renormalized stress-energy 
tensor of the field induced by the detector [8], so that Ai is not observable outside of 
the detector; (c) If Ax was subtracted naively, the uncertainty principle will be violated 
at late times for a is small enough [9]. 

5. Applications to RQI: quantum teleportation 

Quantum teleportation is not only of practical values but also of theoretical interest 
because it contains many illuminating manifestations of quantum physics, clarifying 
fundamental issues such as quantum information and classical information, quantum 
nonlocality and relativistic locality, spacelike correlations and causality, etc. [21 [681 E2] 

The first protocol of quantum teleportation is given by Bennett et al. in [70], where 
an unknown state of a qubit C is teleported from one spatially localized agent Alice to 
another agent Bob using an entangled pair of qubits A and B prepared in one of the 
Bell states and shared by Alice and Rob, respectively. This idea is then adapted to the 
systems with continuous variables by Vaidman [71], who introduces an EPR state [72] 
for the shared entangled pair to teleport an unknown coherent state. Braunstein and 
Kimble (BK) [73] generalized Vaidman's scheme from EPR states with exact correlations 
to squeezed coherent states. In doing so the uncertainty of the measureable quantities 
reduces the degree of entanglement of the y45-pair as well as the fidelity of teleportation. 

To explore how the Unruh effect affects teleportation, Alsing and Milburn made the 
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first attempt of calculating the fidelity of quantum teleportation between two moving 
cavities in relativistic motions [71] - one is at rest (Alice), the other is uniformly 
accelerated (Rob) in the Minkowski frame - though their result is not quite reliable 
[75| I76| 177] . Then Landulfo and Matsas [78] considered quantum teleportation in the 
future asymtotic region in a two-level detector qubit model where Rob's detector is 
uniformly accelerated and interacting with the quantum field only in a finite duration. 
Alternatively, Shiokawa [79] has considered quantum teleportation in the UD detector 
theory with the agents in similar motions but based on the BK scheme in the interaction 
region. More recently, the relativistic effects of quantum information associated with 
the quantum field have been taken into account carefully in [18J, as summarized below. 

Consider a setup with detectors A and C held by Alice, who is at rest in 
space with the worldline z A = z c = (£,1/6,0,0) in the Minkowski frame, and 
detector B held by Rob, who is uniformly accelerated along the worldline z B = 
(a" 1 sinh ar, a -1 cosh ar, 0, 0), < a < 6, where r is Rob's proper time, namely, 
t a = r c = t and r B = r. Suppose the initial state of the combined system at t = r = 
is a product state p§ ®Pab®Pc of the Minkowski vacuum of the field p $x = |0m) (0m| 5 
a two-mode squeezed state p AB of the detectors A and B, and a coherent state of the 
detector C, denoted p^ = \a) c (a\, which is the quantum state to be teleported. To 
concentrate on the best fidelity of quantum teleportation that the entangled AB-p&ir 
can offer, however, we assume the dynamics of p@ is frozen. Also we design p AB so that 
it goes to an EPR state with the correlations while Qa + Qb and Pa — Pb ar< 3 totally 
uncertain as its squeezed parameter r\ — > oo. 

In the Minkowski frame, at t — 0, the detectors A and B start to couple with the 
field, while the detector C is isolated from others. At t = t\ when the reduced state of the 
three detectors continuously evolves to p ABC (K, A; tx), a joint Gaussian measurement 
by Alice is performed locally in space on A and C so that the post-measurement 
state right after t x collapses to p ABC (K,A;^) = p A l(K A , K c , A A , A c )p B (K B , A B ), 
where p" AC {K A 1 K c , A A , A c ) is a two-mode squeezed state of detectors A and C with 
displacement (3 = (5r + iflj, which is the outcome Alice obtains. 

5.1. Entanglement and Pseudo-Fidelities in Different Frames 

Quantum entanglement of the detector pair A and B in a Gaussian state is fully 
determined by the symmetrized two-point correlators Qij, Vij, and IZij, i,j = A,B 
in the quantum state (fT5j) [T6| |80| [HI] . Since the worldlines of the two detectors do not 
intersect, the entanglement between the detectors at some moment will be a kind of 
spacelike correlation and depend on the spatial hypersurface on which the entanglement 
is defined. The time evolution of the entanglement will hence depend on how the 
spacetime is foliated by spacelike hypersurfaces. In general entanglement of the detectors 
is incommensurable with the physical fidelity of quantum teleportation, which is a kind 
of timelike correlation. 

To compare these two correlations, let us first imagine that Rob receives the 
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outcome /3 of Alice's joint measurement and make the proper operation on detector 
B instantaneously at T\(ti) when the worldline of B intersects the ti-slice (see figure [2]). 
According to the outcome (3 obtained by Alice, the operation that Rob should perform 
on detector B is a displacement by (3 in phase space of B from p B to p out , which is 
defined on a time-slice right after the one where the post-measurement state Pabc ls 
defined. The "pseudo-fidelity" of quantum teleportation from \a) c to \a) B is defined 
as F(/3) = B (a \p ou t\ a ) bI^bPouv where Tt B p out = P{0) has been normalized to be the 
probability of finding the outcome /3. Then the averaged pseudo-fidelity is defined by 



F av = 1/2 is known as the best fidelity of "classical" teleportation using coherent states 
[73] , without considering the coupling of the UD detectors with the environment. 

Two results in the ultraweak coupling limit, one in the Minkowski frame, the other 
in the "quasi- Rindler frame" are shown in figure [2] (middle, right). One can see 
that the degree of entanglement (logarithmic negativity Ej^f) evolves smoothly while 
the averaged pseudo-fidelity evaluated in whatever reference frame oscillates in t\ or 
T\. Even at very early times F av drops below 1/2 frequently when E^ is still large. 
Clearly the oscillation of F av here is mainly due to the distortion of the quantum state 
of AS-pair from their initial state (caused by the alternating squeeze-antisqueeze natural 
oscillations of their quantum state) rather than the disentanglement between them. 

In our setup, when t± gets larger, the time dilation of detector B becomes more 
significant and so detector B appears to change extremely slowly in the Minkowski 
frame, while in the quasi-Rindler frame, it is the time-dilation of detector A become 
obvious and so detector A looks frozen for the Rindler observer when T\ gets larger. In 
both cases it makes the frequency of the oscillation of F av approaches Q for larger times. 

The peak values of F av , denoted by F+, fall below the fidelity of classical 
teleportation 1/2 (at the moment denoted ty^ or 7i/ 2 ) always earlier than the dis- 
entanglement time tdE or t^e when Ej^ become zero in both frames. Indeed, in Appendix 
A of [18] it is shown that entanglement between detectors A and B is necessary to provide 
the advantage of quantum teleportation, at least in the ultraweak coupling limit. 

In the Minkowski frame, both F+ v and Eyv are insensitive to a in the ultraweak 
coupling limit. In contrast, the dependence on the proper acceleration a is obvious in the 
quasi-Rindler frame, where the larger a, the earlier T1/2 and the earlier disentanglement 
time TdE- Here we see the frame dependence of entanglement dynamics. 

Beyond the ultraweak coupling limit, both F av and E_\f are strongly affected by 
the environment. In most cases quantum entanglement disappears quickly both in 
the Minkowski frame and the quasi-Rindler frame due to strong interplays with the 
environment, and the averaged pseudo-fidelity F av drops below 1/2 even quicker. 

IT By a quasi-Rindler frame we refer to the coordinate system in which each time-slice almost overlaps 
a Rindler time-slice in the R-wedge but the part in the L-wcdge has been bent to the region with 
positive t to make the whole time-slice located after the initial time-slice for the Minkowski observer, 
as illustrated in figure [2] 




(20) 
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Figure 2. (Left) Setup for quantum teleportation from Alice (thick dotted line) to 
Rob (thick dashed for section 15.11 and thick solid for I5.2[) . The gray solid curve and 
the horizontal line represent a r-f-slice in the quasi-Rindler frame and a ti-slice in the 
Minkowski frame, respectively. The shaded region is in the future lightcone of the 
measurement event by Alice, and the hypersurface t = x 1 is the event horizon of Rob 
for r 2 — > oo. (Middle) A comparison of the averaged pseudo-fidelity F av (solid) in the 
Minkowski frame, the correlator (Q 2 _)/20 (dotted), Q_ = Qa~Qb and the logarithmic 
negativity Ej^/3.5 (dot-dashed) at early times in the ultraweak coupling limit. (Right) 
F av in the quasi-Rindler frame with the same parameters (solid) and with different a 
(dotted, gray). The upper and middle dashed lines represent F+, and 1/2. 

5.2. Physical fidelity and "entanglement on the lightcone" 

Suppose Rob stops accelerating at his proper time T2, after this moment Rob moves 
with constant velocity, while Alice stays at rest and performs the joint measurement on 
A and C at t\ (see figure EJ). In this setup the classical information from Alice traveling 
at the speed of light can always reach Rob, though the acceleration of detector B is not 
uniform - For the dynamics of the correlators similar situations we refer to Ref. |14j . 
more results can be found in |82j . 

Suppose Rob performs the local operation on B at some moment t p > r^ dv when he 
received the information traveling in lightspeed from Alice (see figure [2]for the definition 
of rf ). Then, similar to (120]) . the averaged physical fidelity should be given by 
Fav = J d 2 (3 B (a \ Poutijp) \a) B , where p out (T P ) is obtained by performing a displacement 
on P b (t p ) which started with the post-measurement state P b {t\) with n = a~ l sinh -1 at\ 
(when the quantum state collapses in the Minkowski frame) and evolves from Ti to r p . 
Nevertheless, an analysis similar to p2] shows that the correlators in the reduced state 
of detector B observed in all reference frames will become the same collapsed ones at 
the moment when detector B is entering the future lightcone of the measurement event 
by Alice [r B = rf ). So we are allowed to collapse the wave functional on a time-slice 
intersecting Alice's worldline at r A = t\ and Rob's at r B = r"* — e, e — > 0+. 

If we further assume that mutual influences are small and Rob performs the local 
operation at t p = rf +e right after the classical information from Alice is received, then 
the continuous evolution of the reduced state of detector B from r M to r P is negligible, 
and so we can directly compare the physical fidelity at r^ dv + e with entanglement 
between detectors A at t\ and B at r^ dv — e. Again from Appendix A of [18] with the 
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Figure 3. F+, — 1/2 (solid) and Ej^/\Q (dashed) in section 15721 as a lunction of the 
moment of the joint measurement t\. The black curves are those physical fidelities and 
the degrees of entanglement evaluated "on the lightcone" . while the gray curves are 
those pseudo-fidelities and the degrees of entanglement evaluated on the ti-slice in the 
Minkowski frame. 



proper time of detector B substituted by r^ dv (actually r^ dv ± e, e — > 0+), quantum 
entanglement of AB-p&ir evaluated almost on the future lightcone of the measurement 
event by Alice is still a necessary condition of the best averaged physical fidelity of 
quantum teleportation beating the classical one in the ultraweak coupling limit. 

The number of peaks of the physical F av in the same duration of t\ in this more 
realistic case is much more than the one for the averaged pseudo-fidelity, because it takes 
a long time from Ti(t\) to the moment rf when the classical signal from Alice reaches 
Rob, during which detector B has oscillated for many times. In figure [3] we see that 
the moment t\ = ti/2 when the best averaged physical fidelity of quantum teleportation 
F+ v drops to 1/2 is earlier than any F+ v of pseudo-fidelity has. The larger ar 2 , the later 
T adv -p^k an( j so lower value of the physical F+ v at that time due to the longer 
time of coupling with environment. When aT 2 is large enough, Tf dv is so large that ti/ 2 
is almost the moment that Alice enters the event horizon of Rob for r 2 — > oo. 

6. Outlook 

Detector-field interaction is a very useful and versatile system to expound and explore 
many known and unknown effects of relativistic quantum information. Continuing the 
vein of quantum teleportation as example, it is not difficult to generalize the systematics 
to curved spacetimes, from weak gravitational field as in the Earth's environment to 
black hole spacetimes. 

In another vein, learning from the new techniques and ideas in the study of detector- 
field interaction as was done for the Unruh effect [HI E] , and using the well-known corre- 
spondence between the Rindler and the Schwarzschild spacetimes, one may go beyond 
the test-field description of black hole physics and study how backreaction from the field 
impacts on the evolution of the black hole and the "information loss" issues. 
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